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Abstract 

A system of two particles with spin s = and s = 77 respectively, moving in a plane 

is considered. It is shown that such a system with a nontrivial spin-orbit interaction can 

""^ ■ allow an 8 dimensional Lie algebra of first-order integrals of motion. The Pauli equation 

is solved in this superintegrable case and reduced to a system of ordinary differential 

JL( equations when only one first-order integral exists. 
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I INTRODUCTION 

A Hamiltonian system with n degrees of freedom in classical mechanics is called integrable if 
it allows n functionally independent integrals of motion {Xi, . . . , X n }. These integrals must be 
well-defined functions on phase space and be in involution. The Hamiltonian H belongs to this 
set of n integrals. A "superintegrable system" is one that allows some additional integrals of 
motion, {Y 1; . . . , Y k } such that the set {Xi, . . . , X n , Yi, . . . , Y k } is functionally independent. The 
integrals {Yi, . . . , Y k } are not necessarily in involution among each other, nor with the X^. A 
system is maximally superintegrable if we have k = n — 1, minimally superintegrable for k — 1. 

The concepts of integrability and superintegrability are also introduced in quantum me- 
chanics. The only difference is that the integrals of motion are now well-defined linear quantum 
mechanical operators, assumed to be algebraically independent. 1,6_9 ' 12 ' 13 ' 15 ' 16 ' 18 ' 19 ' 21 

The best known superintegrable systems are the Kepler, or Coulomb system 1 ' 8 and the 
harmonic oscillator. 12,16 They are characterized by the fact that all finite classical trajectories 
in these systems are periodic. In quantum mechanics these systems are exactly solvable, i.e. 
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their bound state energy levels can be calculated algebraically and their wave function expressed 
in terms of polynomials. 

The above properties are shared by all other known maximally superintegrable systems (see 
e.g. 19 ). 

A systematic search for superintegrable systems and their properties was started quite some 
time ago. 6 ' 7 ' 9 ' 15,21 Originally the approach concentrated on Hamiltonians of the type 

H = - 1 -A + V(r) (1) 

in 2- and 3- dimensional Euclidean spaces with the restriction that all integrals of motion should 
be first- or second-order polynomials in the momenta. More recently the study of superintegrable 
systems with second-order integrals of motion was extended to curved spaces and also higher- 
dimensional ones (see 13 for some recent results and an extensive list of references). 

For Hamiltonians of the type ([1]) with second-order integrals of motion there is a close relation 
between integrability and the separation of variables in the Schrodinger and Hamilton- Jacobi 
equations. Typically superintegrable systems of this type are multiseparable: they allow the 
separation of variables in more than one system of coordinates. 

This relationship between integrability and separability breaks down in other cases. Thus 
for natural Hamiltonians of the type (CQ) the existence of third-order integrals of motion does 
not lead to the separation of variables. 4 ' 5 ' 10 ' 11 Furthermore, if we consider velocity dependent 
potentials (e.g. related to magnetic field), 

H=-^A + V(r) + (A,p), (2) 

then quadratic integrability no longer implies the separation of variables. 2, 3,17 ' 2 ° 

In this article we initiate the study of integrability and superintegrability in a different 
type of system, namely one involving particles with spin. More specifically, we consider two 
nonrelativistic quantum particles, moving in a plane, one with spin |, the other with spin 0. In 
this case the Hamiltonian will be a matrix operator, acting on two component spinors and we 
can decompose it in terms of the identity matrix / and the Pauli matrices <j{ (i = 1, 2, 3). 
From the physical point of view the most interesting Hamiltonian to consider would be 



2m 



A + Vo(rO + £{fi(rO, (*,£)} (3) 



in the 3- dimensional Euclidean space £3 (we drop the matrix / whenever this does not cause 
confusion). The curly bracket represents an anticommutator and V a (f), a = 0, 1 are real func- 
tions. The Hamiltonian is Hermitian and satisfies the requirements of parity and time reversal 
invariance. The spin-orbital interaction term V\_(r)(a,L) is the standard one in quantum me- 
chanics. 14 

In this paper, the first one on integrability and superintegrability for particles with nonzero 
spin, we restrict to the 2- dimensional space E 2 . The orbital angular momentum L then only 
has one component, L 3 , perpendicular to the xy- plane E 2 . The scalar product (a, L) reduces 
to 03L3 (since L\ and L 2 are zero). We shall set the reduced mass m of the two particle system 
equal to m = 1 and use units in which the Planck constant is % = 1 (we do not need to consider 
a classical limit here). Finally, the Hamiltonian to be considered in this article is: 

H = ^(Pi 2 +P2 2 ) + V {x,y) + V 1 {x,y)a 3 L 3 + -a^L^V^x.y)) (4) 



with 

Pi = -id x , p 2 = -id y , L 3 = i(yd x - xd y ), a 3 = I " " 1 . (5) 




A priori the functions Vo(x, y) and Vi(x,y) are arbitrary. In specific physical applications 
they may be related or they may both be specified. Our aim is to determine the conditions on 
these two functions, under which one or more integrals of motion exist. 

We request that at least one first-order integral of motion should exist, namely 

X = (A oPl + B oP2 + O )I + (A lPl + B lP2 + <t> x )(T3 , (6) 

where A^, B^ and M (/i = 0, 1) are real functions of x and y. These functions as well as the 
potentials Vo and V\ are to be determined from the commutativity condition 

[H, X} = 0. (7) 

The general formulation is set up in Section [TT1 where we determine A^ and B^ and obtain 
the partial differential equations (PDE) that (f>o,(j)i,Vo and Vi must satisfy. In Section HTT1 we 
consider a special case when the Hamiltonian (pD) allows 6 independent nontrivial integrals of the 
type ([6]). They generate an 8-dimensional symmetry group of the system. Section HVl is devoted 
to more general integrable Hamiltonians, allowing just 1 first-order integral. The system of 
equations 

H^ = EV, Xm = A^ (8) 

is studied in Section [V] We shall call the system ([8]) the Pauli system. Some conclusions and 
future directions are outlined in the final Section IVII 

II FORMULATION OF THE PROBLEM 

In order to obtain determining equations for the coefficients A^, B^ and M , (/i = 0, 1) in the 
integral ([6]) we impose the commutation relation (171). The commutator will involve terms of 
the type pf, p 2 a 3 , p 2 2 , p 2 2 a 3 , p x p 2 , Pip 2 <r 3 , Pi, Pi&3, Pi, P2O3, I and a 3 . We have a 2 = I, so no 
higher powers of a 3 appear. We set the coefficients of each of the above terms equal to zero. 
This gives us 12 linear partial differential equations for A^, B^ and M . Those coming from the 
coefficients of second-order terms in the momentum imply that A^ and B^ are linear functions 
and we obtain, for any potentials Vo and V\ 

Ap, = u,jy + a^, B^ = -uj^x + bf,, (9) 

where u^, a^ and b^ are real constants. The coefficients of p, pa 3 , I and a 3 in the commutator 
provide an overdetermined system of six first-order PDEs for the four functions Vo, Vi, 0o an d 
0i. They are 

M ,x = ^,i-v[-b v Vi - {u v y + a v )yV hx + {u v x - K)yV liy ] , 

(10) 
M ,y = 5»,i-u[a v Vi + (w v y + a v )zV 1)X - (v„x - b u )xV lt y] , 



(u^y + ajVoj. + {-UpX + b^)V 0:y = 5^ v {x(f)^ y - y<j> v ^)V x , Qu, v = 0, 1) . (11) 

The coefficients of I and 03 a priori involve second-order derivatives of V±(x,y). These 
second-order terms cancel, once equations © and (TTOT) are taken into account. This leads to 
the two first-order equations ffTTl) . 

Before solving this system, let us introduce "allowed transformations" that leave the Hamil- 
tonian (j3j) form invariant, i.e. change only the functions Vo(x, y) and Vi(x,y). Such transfor- 
mations will be used to simplify Hamiltonians, integrals of motion and also the equations to be 
solved. 

Allowed transformations for any potentials Vq and V% are: 

1. Rotations in the xy-plane. 

2. Gauge transformations of the form 

H=U~ 1 HU, U=( 6 ™ e ° iQ V a = «(£), £=|. (12) 

The transformation of the potentials is 

A, y = V4 + (l + 4)(i^ 



^1 = ^ + ^, V = V + (l + ^)(-— + aV 1 ). (13) 



For certain specific potentials V\ further allowed transformations exist (for any Vo), namely 
simultaneous translations and gauge transformations 

1. V\ = 7 = const 

The allowed transformations are given by 

x = x + x , y = y + y , a = -y(y x - x y) (14) 

and the transformation of the potentials is 

Vi = Vi=7, V (x,y) = V {x + x ,y + y ) - -j 2 (x 2 + y 2 + 2(xx + yy )) . (15) 

2. Vi = Vi(x) 

The allowed transformation is 

V 1 = V 1 , V (x, y) = V (x, y + y ) - ^i 2 y (?/o + 2y) . (16) 



Let us now return to equation (TlOl) . The compatibility conditions <j)^ xy = 4>^,yx imply 

(u^y + a^)(xV hxx + yV hxy + 3V ljX ) + (-u^x + b^yV^yy + xVi,^ + 3Vi j2/ ) = , (/i = 0, 1) .(17) 

In general, ( iT7j) represents an overdetermined system of two different equations for the 
potentials Vi(x, y) and this system can be written as 

xV 1>xx + yV 1>xy + 3V ljX = , 
yV lM + xV 1<xy + 3V 1>y = . (18) 



An exception occurs if the two equations (1171) coincide. This happens if the constants figuring 
in equation (flTl) satisfy the three following equations: 

uj Q ai — ujiclq = 0, u>obi — cuibo = 0, aib — a bi = . (19) 

We shall treat the case ( ITBl in Section [III] below. The case when (TT9l) is satisfied and Vi(x, y) 
satisfies only one equation ( TPTj) will be considered in Section [TV] 

III SPIN ORBITAL INTERACTION WITH 
KINEMATICAL INVARIANCE GROUP 

Let us now solve equations ( TT8l) . We transform the first equation to characteristic variables, 
solve and substitute into the second equation. The result is 

V 1 (x,y) = y+^, £ = £, y = const . (20) 

Jb Jb 

Comparing with equation f|T3|) we see that we can annul the function G(£) by a gauge transfor- 
mation. Thus we have 

V 1 = y. (21) 

Substituting ( 12T1) into equations ( TTUj) and ( TTT1) we obtain 

K) = 2^ 2 ( x2 + ?/ 2 )' 0o = ~7(M-ai1/), 0i = -j(b x - a y) . (22) 

The Hamiltonian thus has the form 

# = -^A + l 7 V + ?/ 2 ) + 7^3, l^R- (23) 

Since H does not depend on the constants cj m , a M and b^ we obtain 6 independent integrals of 
motion, generating the symmetry group of this Hamiltonian. 
A basis for the symmetry algebra is given by the 8 operators 



L± ~- 


= i(yd x -xd y )(I±a 3 


x± -- 


= (id x ^yy)(I±a 3 ), 


Y± ~- 


= (id y ± 7x)(J ± cr 3 ) , 


I± ~- 


= I±a 3 . 



(24) 
The nonzero commutation relations are 

[L ± , X ± ] = 2iY ±) [L ± , Kt] = -2iXt, [X ± , Y ± ] = ±4z 7 J ± . (25) 

The symmetry algebra is thus isomorphic to the direct sum of two central extensions of the 
Euclidean Lie algebra 

L ~ e+(2) © e_(2) = {L + , A + , F + , /+} © {L_,X_, F_, J_} . (26) 

The Casimir operators of L are 

C± = X ± 2 + Y ± 2 ± 4 7 L ± /± , 7± = / ± 0-3 (27) 



and we have 



H=\(C + + C^). (28) 



The integral of motion X is a linear combination of the 8 operators (|24"|) with arbitrary real 
constant coefficients. Such operators X can be classified into conjugacy classes under the action 
of the group generated by the algebra (J241) . The conjugacy classes that lead to different types 
of solutions of the Pauli system (jSJ) can be represented by 

Xi = L + + aL_, X 2 = L + + al_, X 3 = X + + aX_ , a e i? . (29) 

The Hamiltonian (j2"3"l) is not only integrable, but actually "first-order superintegrable" . For 
particles of spin first-order superintegrability occurs only for free motion. Notice that if we 
set the spin-orbit interaction equal to zero in equation (T231 (i.e. 7 = 0), we obtain free motion. 

IV HAMILTONIANS ALLOWING ONE 
FIRST-ORDER INTEGRAL 

Let us now consider the case when equations (fT9l) are satisfied. The two equations (IT7j) then 
coincide and the potential Vi(x,y) satisfies just one second-order PDE. The equation (ITT)) is of 
hyperbolic type. Its characteristic variables are 

£=-, r) = -uii{x 2 + y 2 ) -bix + aiy. (30) 

x 2 

Here we shall just consider two interesting special cases. 

a) u>i 7^ 0, a = b = 0, d\ = b\ = 

We transform equation (fT7|) to polar coordinates and obtain 



pVi, P e + 2Vi,e = , x = pcos6, y = psm6 . (31) 

The potential Vi hence is 

^1 = f(p) + \g(0) , (32) 

p 

where f(p) and g(8) are arbitrary. Comparing with equation ffTBT) we see that the function 
p(0) can be set equal to g = by a gauge transformation. Solving ( 1101) and (1111) we obtain 

V = V (p), V, = V 1 (p), O = 1 = O, (33) 

X = {ujo + uj^L^. (34) 

b) u = ui = 0, a x 2 + 61 2 ^ 

Equation (|T7|) in characteristic variables (1311 is 

K,fc + ~V U = (35) 



and we obtain 

V, = Fi(oij/ - hx) + ^— . (36) 

(aiy - hx) 2 

By a gauge transformation we set F 2 = and rotate in the xy- plane to obtain 

V x = V x {x). (37) 

From equations ffTUl) and ffTTl) we obtain 

K, = y Vi 2 (x) + F(x) , (38) 

(f) = -b x [Vi(x)dx, (p! = -b jVi{x)dx. (39) 



Let us put b\ = 1, bo ^ 0. We then have 

0i = &o</>o , V x {x) = -(po(x) , 

Vo(x,y) = ^y 2 [(f> (x)} 2 + F(x) (40) 

and 

X = -ib d y + (p (x) + {-id y + b (po(x))cr 3 . (41) 

We shall call the case (1331) the "polar" case, (137|) . (|38|) the "Cartesian" one, because of the form 
of the operator X in f |34|) and f j4Tl) . respectively. 

V SOLUTIONS OF THE PAULI EQUATION 

In this section we shall analyze and solve the pair of equations ([8]) for the different superinte- 
grable, or first-order integrable cases, found above. 

1. The superintegrable Hamiltonian 

Let us consider the Hamiltonian ( |23l) with 7^0, i.e. a constant spin-orbital potential and 
a harmonic oscillator spin-independent one. The Hamiltonian commutes with the entire kine- 
matical algebra (1241) . We shall choose the operator X of equation ([8]) in the form of one of the 
different one-dimensional subalgebras shown in equation (129]) and consider each of the three 
cases separately. The potentials V (x,y) and Vi(x,y) in this case have the form (|33j) and (T40|) 
simultaneously. Hence we can separate variables in polar coordinates, as well as in Cartesian 
ones. Moreover, we can consider a mixed case: separation in different coordinate systems for 
the upper and lower components of ^/. 

a) Polar case 

We introduce polar coordinates (p, 9) and choose the operator X in the form 

X = ^(o a 2 ) de > a ^°- (42) 



The condition X^f = A\& provides a wave function in the form 

n^)=( F ; { l p) fil)- («) 

\ F 2 {p)e °2 J 

Substituting into the Pauli equation with Hamiltonian ( |23l) we find that the function Fi(p) 
satisfy 

f l f d 2 Id 1 X\ l n ] r _ A.^ 
[ 2 v <9p pc*p p l ai l > 2 J a» 

so both components satisfy radial harmonic oscillator type equations. 
The solution of equation (144")) is 

Fi(p) = JV»™ e-^ 2 p^L^( 1P 2 ) , (45) 

where /^'(z) are Laguerre polynomials. The quantum number A satisfies 

A A 

— = m 1 , — —m 2 , (46) 

ai a 2 



hence we must choose 



rational. The energy satisfies 



a 2 _ m 1 
a% m 2 



(47) 



E — 7ml = j(2ni + \mi\ + 1) , 

E + -?m 2 = ~f(2n 2 + \m 2 \ + 1) , (48) 

so the two radial quantum numbers are constrained by 

2(n 2 - m) = mi f— + l) - \mx\ (\—\ - l) • (49) 

\a 2 J \ a 2 J 

If we normalize to have 

/ / (\V 1 \ 2 + \* 2 \ 2 )pdpd9=l, (50) 

Jo Jo 

we must put the normalization constants in ( 145]) equal to 



Ly\mi\+1 / ^J 

V 27T V (n» + |77li|)! 



b) Cartesian case 

We choose the operator to be diagonalized in the form 



-H ai(! V 7i) a2 K - 7 ,)). **°- ^ 



The equation X^ = A\l/ implies 



F 1 (x)e-^ {X - ailx)y 



* ! ^ : _^-.,,-.,,„ 1 • ««' 



F 2 (x)e-^ (A+ " 2 "''" 

Substituting into the Pauli equation if \1/ = i?\l/ with H as in equation (123]) we obtain 

A ^ 2 



F,-4 7 Mx T - F l + 2EF l = 0, i = 1, 2 (54) 

V 2^7/ 



and hence 



F i = N tH e-^H n .(^x i 



E = 27(ni + -) = 27(n 2 + -) , 7^ = n 2 = n , 



AT n 



Xi, 2 = X =F - , (55) 



where H n (z) is a Hermite polynomial. 

c) Mixed case 

Let us take the operator X in the form 



The wave function will then be 

/ F x {p)e^ e \ 

^ = \F 2 {x)e-^ X+a ^) ' (57) 

where (p, 6) are polar coordinates, (x, y) Cartesian ones. The function F 1 (p) will satisfy equation 
(|44|) with i = 1, F 2 (x) equation (154)) with z = 2. We hence obtain 

F 2 (x) = N n2 e~^H n2 (^x 2 ), (58) 

with 

£ = 7(7711 + 2ni + |toi| + 1) = 2 7 (n 2 + -) . (59) 



2. Hamiltonians with one first order integral 

a) Polar case 

We consider the potential Vq = Vo(p), V\ = V\(p) as in equation ( J33l . We write the integral 
( |34l) in the form 

X = -i( ^ ^ \ d e , a, ^ (60) 

and the equation X\l/ = A\l/ implies 

*te»)=|'^ e !t|. («) 

F 2 {p)e a 2 

Substituting into the Pauli equation, we find that the radial functions iq, F2 satisfy: 

HS + ^-7^) + ^ ± ^} F ' = ^' < = 1 ' 2 - (62) 

For instance, choosing 

V = -, V 1 = 4 , 63 

P P 

we can solve equation (162]) in terms of Coulomb wave functions. 

b) Cartesian case 

Let us now consider Vo, V\ and the integral X as in (1401) . (141]) . The equation X\l/ = A\I/ with X 
as in (J4"I]) implies 

V X2W e "o-i y 

The Pauli equation then reduces to the following two ODEs: 

.. f (A T (fr o ±l)0 o ) 2 ^ ^, At fa ±l)fo , orv v 9 J n • 19 (M 

Xi ~) (b ±l) 2 ± 2x( ^o , h ZjTT/j + 2F(x) -2£^x; = 0, 2 = 1,2. (65) 

To solve these equations, or analyze further, we would have to specify the two functions F(x), 

(f) (x). 

VI CONCLUSIONS 

Let us first of all compare the problem of integrability and superintegrability for particles with 
spin s = and spin s 7^ 0, in this case s = ~. The spin zero case with a scalar potential in the 
two-dimensional Euclidean space Ei corresponds to the Hamiltonian (CQ). First order superinte- 
grability is trivial: it requires V(x) = const, and corresponds to free motion. Superintegrability 

10 



with one first-order and one second-order integral occurs for the potentials V = ar 2 , ar _1 , ax 
and ax~~ 2 . 18 Quadratic superintegrability leads to 4 families of potentials, 9 ' 21 each depending 
on 3 significant constants and allowing the separation of variables in at least two coordinate 
systems. 

By contrast for s = | first-order superintegrability leads to a nontrivial system, namely 
the Hamiltonian (1231) with the symmetry algebra (124")) . The Hamiltonian allows the separation 
of variables in polar coordinates (see (j!5|) ) and "R-separation" in Cartesian ones (see (1531 ). 
Indeed, in equation (153]) there is a term involving the product xy, that does not depend on the 
separation constant A. The same is true for Hamiltonians allowing just one first-order integral. 
In the polar case we have separation (see (I6T1) ). in the Cartesian one R-separation ( 1641) . 

The next step in this research program will be to look for integrable and superintegrable 
systems with spin in Euclidean 3-space. This would provide a realistic and solvable model for 
pion-nucleon and possibly nucleon-nucleon interactions. 
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